Variants in Human ATP Synthase Mitochondrial Genes: Biochemical Dysfunctions, Associated Diseases, and Therapies.

Mitochondrial ATP synthase (Complex V) catalyzes the last step of oxidative phosphorylation and provides most of the energy (ATP) required by human cells. The mitochondrial genes MT-ATP6 and MT-ATP8 encode two subunits of the multi-subunit Complex V. Since the discovery of the first MT-ATP6 variant in the year 1990 as the cause of Neuropathy, Ataxia, and Retinitis Pigmentosa (NARP) syndrome, a large and continuously increasing number of inborn variants in the MT-ATP6 and MT-ATP8 genes have been identified as pathogenic. Variants in these genes correlate with various clinical phenotypes, which include several neurodegenerative and multisystemic disorders. In the present review, we report the pathogenic variants in mitochondrial ATP synthase genes and highlight the molecular mechanisms underlying ATP synthase deficiency that promote biochemical dysfunctions. We discuss the possible structural changes induced by the most common variants found in patients by considering the recent cryo-electron microscopy structure of human ATP synthase. Finally, we provide the state-of-the-art of all therapeutic proposals reported in the literature, including drug interventions targeting mitochondrial dysfunctions, allotopic gene expression- and nuclease-based strategies, and discuss their potential translation into clinical trials.

Structural and biochemical characterization of Arabidopsis alcohol dehydrogenases reveals distinct functional properties but similar redox sensitivity.

Alcohol dehydrogenases (ADHs) are a group of zinc-binding enzymes belonging to the medium-length dehydrogenase/reductase (MDR) protein superfamily. In plants, these enzymes fulfill important functions involving the reduction of toxic aldehydes to the corresponding alcohols (as well as catalyzing the reverse reaction, i.e., alcohol oxidation; ADH1) and the reduction of nitrosoglutathione (GSNO; ADH2/GSNOR). We investigated and compared the structural and biochemical properties of ADH1 and GSNOR from Arabidopsis thaliana. We expressed and purified ADH1 and GSNOR and determined two new structures, NADH-ADH1 and apo-GSNOR, thus completing the structural landscape of Arabidopsis ADHs in both apo- and holo-forms. A structural comparison of these Arabidopsis ADHs revealed a high sequence conservation (59% identity) and a similar fold. In contrast, a striking dissimilarity was observed in the catalytic cavity supporting substrate specificity and accommodation. Consistently, ADH1 and GSNOR showed strict specificity for their substrates (ethanol and GSNO, respectively), although both enzymes had the ability to oxidize long-chain alcohols, with ADH1 performing better than GSNOR. Both enzymes contain a high number of cysteines (12 and 15 out of 379 residues for ADH1 and GSNOR, respectively) and showed a significant and similar responsivity to thiol-oxidizing agents, indicating that redox modifications may constitute a mechanism for controlling enzyme activity under both optimal growth and stress conditions.

Genetic variants affecting NQO1 protein levels impact the efficacy of idebenone treatment in Leber hereditary optic neuropathy.

Idebenone, the only approved treatment for Leber hereditary optic neuropathy (LHON), promotes recovery of visual function in up to 50% of patients, but we can neither predict nor understand the non-responders. Idebenone is reduced by the cytosolic NAD(P)H oxidoreductase I (NQO1) and directly shuttles electrons to respiratory complex III, bypassing complex I affected in LHON. We show here that two polymorphic variants drastically reduce NQO1 protein levels when homozygous or compound heterozygous. This hampers idebenone reduction. In its oxidized form, idebenone inhibits complex I, decreasing respiratory function in cells. By retrospectively analyzing a large cohort of idebenone-treated LHON patients, classified by their response to therapy, we show that patients with homozygous or compound heterozygous NQO1 variants have the poorest therapy response, particularly if carrying the m.3460G>A/MT-ND1 LHON mutation. These results suggest consideration of patient NQO1 genotype and mitochondrial DNA mutation in the context of idebenone therapy.

Determination and Kinetic Characterization of a New Potential Inhibitor for AmsI Protein Tyrosine Phosphatase from the Apple Pathogen Erwinia amylovora.

Erwinia amylovora is a Gram-negative bacterium, responsible for the fire blight disease in Rosaceae plants. Its virulence is correlated with the production of an exopolysaccharide (EPS) called amylovoran, which protects the bacterium from the surrounding environment and helps its diffusion inside the host. Amylovoran biosynthesis relies on the expression of twelve genes clustered in the ams operon. One of these genes, amsI, encodes for a Low Molecular Weight Protein Tyrosine Phosphatase (LMW-PTP) called EaAmsI, which plays a key role in the regulation of the EPS production pathway. For this reason, EaAmsI was chosen in this work as a target for the development of new antibacterial agents against E. amylovora. To achieve this aim, a set of programs (DOCK6, OpenEye FRED) was selected to perform a virtual screening using a database of ca. 700 molecules. The six best-scoring compounds identified were tested in in vitro assays. A complete inhibition kinetic characterization carried out on the most promising molecule (n-Heptyl ß-D-glucopyranoside, N7G) showed an inhibition constant of 7.8 ± 0.6 µM. This study represents an initial step towards the development of new EaAmsI inhibitors able to act as antibacterial agents against E. amylovora infections.

Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

Design, Pharmacological Characterization, and Molecular Docking of Minimalist Peptidomimetic Antagonists of α4ß1 Integrin.

Integrin receptors mediate cell-cell interactions via the recognition of cell-adhesion glycoproteins, as well as via the interactions of cells with proteins of the extracellular matrix, and upon activation they transduce signals bi-directionally across the cell membrane. In the case of injury, infection, or inflammation, integrins of ß2 and α4 families participate in the recruitment of leukocytes, a multi-step process initiated by the capturing of rolling leukocytes and terminated by their extravasation. In particular, α4ß1 integrin is deeply involved in leukocyte firm adhesion preceding extravasation. Besides its well-known role in inflammatory diseases, α4ß1 integrin is also involved in cancer, being expressed in various tumors and showing an important role in cancer formation and spreading. Hence, targeting this integrin represents an opportunity for the treatment of inflammatory disorders, some autoimmune diseases, and cancer. In this context, taking inspiration from the recognition motives of α4ß1 integrin with its natural ligands FN and VCAM-1, we designed minimalist α/ß hybrid peptide ligands, with our approach being associated with a retro strategy. These modifications are expected to improve the compounds' stability and bioavailability. As it turned out, some of the ligands were found to be antagonists, being able to inhibit the adhesion of integrin-expressing cells to plates coated with the natural ligands without inducing any conformational switch and any activation of intracellular signaling pathways. An original model structure of the receptor was generated using protein-protein docking to evaluate the bioactive conformations of the antagonists via molecular docking. Since the experimental structure of α4ß1 integrin is still unknown, the simulations might also shed light on the interactions between the receptor and its native protein ligands.

The structure of the high-affinity nickel-binding site in the Ni,Zn-HypA•UreE2 complex.

The maturation pathway for the nickel-dependent enzyme urease utilizes the protein UreE as a metallochaperone to supply Ni(II) ions. In Helicobacter pylori urease maturation also requires HypA and HypB, accessory proteins that are commonly associated with hydrogenase maturation. Herein we report on the characterization of a protein complex formed between HypA and the UreE2 dimer. Nuclear magnetic resonance (NMR) coupled with molecular modelling show that the protein complex apo, Zn-HypA•UreE2, forms between the rigorously conserved Met-His-Glu (MHE motif) Ni-binding N-terminal sequence of HypA and the two conserved His102A and His102B located at the dimer interface of UreE2. This complex forms in the absence of Ni(II) and is supported by extensive protein contacts that include the use of the C-terminal sequences of UreE2 to form additional strands of ß-sheet with the Ni-binding domain of HypA. The Ni-binding properties of apo, Zn-HypA•UreE2 and the component proteins were investigated by isothermal titration calorimetry using a global fitting strategy that included all of the relevant equilibria, and show that the Ni,Zn-HypA•UreE2 complex contains a single Ni(II)-binding site with a sub-nanomolar KD. The structural features of this novel Ni(II) site were elucidated using proteins produced with specifically deuterated amino acids, protein point mutations, and the analyses of X-ray absorption spectroscopy, hyperfine shifted NMR features, as well as molecular modeling coupled with quantum-mechanical calculations. The results show that the complex contains a six-coordinate, high-spin Ni(II) site with ligands provided by both component proteins.

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain-domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.

NDUFS3 knockout cancer cells and molecular docking reveal specificity and mode of action of anti-cancer respiratory complex I inhibitors.

Inhibition of respiratory complex I (CI) is becoming a promising anti-cancer strategy, encouraging the design and the use of inhibitors, whose mechanism of action, efficacy and specificity remain elusive. As CI is a central player of cellular bioenergetics, a finely tuned dosing of targeting drugs is required to avoid side effects. We compared the specificity and mode of action of CI inhibitors metformin, BAY 87-2243 and EVP 4593 using cancer cell models devoid of CI. Here we show that both BAY 87-2243 and EVP 4593 were selective, while the antiproliferative effects of metformin were considerably independent from CI inhibition. Molecular docking predictions indicated that the high efficiency of BAY 87-2243 and EVP 4593 may derive from the tight network of bonds in the quinone binding pocket, although in different sites. Most of the amino acids involved in such interactions are conserved across species and only rarely found mutated in human. Our data make a case for caution when referring to metformin as a CI-targeting compound, and highlight the need for dosage optimization and careful evaluation of molecular interactions between inhibitors and the holoenzyme.

PacDOCK: A Web Server for Positional Distance-Based and Interaction-Based Analysis of Docking Results.

Molecular docking is a key method for structure-based drug design used to predict the conformations assumed by small drug-like ligands when bound to their target. However, the evaluation of molecular docking studies can be hampered by the lack of a free and easy to use platform for the complete analysis of results obtained by the principal docking programs. To this aim, we developed PacDOCK, a freely available and user-friendly web server that comprises a collection of tools for positional distance-based and interaction-based analysis of docking results, which can be provided in several file formats. PacDOCK allows a complete analysis of molecular docking results through root mean square deviation (RMSD) calculation, molecular visualization, and cluster analysis of docked poses. The RMSD calculation compares docked structures with a reference structure, also when atoms are randomly labelled, and their conformational and positional differences can be visualised. In addition, it is possible to visualise a ligand into the target binding pocket and investigate the key receptor-ligand interactions. Moreover, PacDOCK enables the clustering of docking results by identifying a restrained number of clusters from many docked poses. We believe that PacDOCK will contribute to facilitating the analysis of docking results to improve the efficiency of computer-aided drug design.

Relevance of ARID1A Mutations in Endometrial Carcinomas.

Since the Cancer Genome Atlas (TCGA) project identified four distinct groups based on molecular alterations, mutation analyses have been integrated into the characterization of endometrial carcinomas (ECs). ARID1A seems to be the subunit more involved in the loss of function of the SWI/SNF complex in ECs. The aim of this study is to define the relevance of ARID1A alterations in a cohort of EC, studying the possible associations between DNA mutation (genomic level), RNA expression (transcriptomic level), and protein expression (proteomic level). A total of 50 endometrial carcinomas were characterized for ARID1A mutations (using targeted DNA next-generation sequencing-NGS), ARID1A gene expression (using RNAseq and qRT-PCR), and ARID1A protein expression (using immunohistochemistry-IHC). Moreover, we have investigated if ARID1A mutations may alter the protein structure, using the Protein Data Bank sequence. We found a good correlation between ARID1A mutations and protein immunostaining, even if we did not find statistically significant differences in the ARID1A expression levels. In conclusion, our data demonstrated that the molecular characterization of ARID1A should be associated with IHC analysis, mainly in those cases harboring "novel" ARID1A mutations or in those alterations with "uncertain" pathogenic significance.

Cloning the barley nec3 disease lesion mimic mutant using complementation by sequencing.

Disease lesion mimic (DLM) or necrotic mutants display necrotic lesions in the absence of pathogen infections. They can show improved resistance to some pathogens and their molecular dissection can contribute to revealing components of plant defense pathways. Although forward-genetics strategies to find genes causal to mutant phenotypes are available in crops, these strategies require the production of experimental cross populations, mutagenesis, or gene editing and are time- and resource-consuming or may have to deal with regulated plant materials. In this study, we described a collection of 34 DLM mutants in barley (Hordeum vulgare L.) and applied a novel method called complementation by sequencing (CBS), which enables the identification of the gene responsible for a mutant phenotype given the availability of two or more chemically mutagenized individuals showing the same phenotype. Complementation by sequencing relies on the feasibility to obtain all induced mutations present in chemical mutants and on the low probability that different individuals share the same mutated genes. By CBS, we identified a cytochrome P450 CYP71P1 gene as responsible for orange blotch DLM mutants, including the historical barley nec3 locus. By comparative phylogenetic analysis we showed that CYP71P1 gene family emerged early in angiosperm evolution but has been recurrently lost in some lineages including Arabidopsis thaliana (L.) Heynh. Complementation by sequencing is a straightforward cost-effective approach to clone genes controlling phenotypes in a chemically mutagenized collection. The TILLMore (TM) collection will be instrumental for understanding the molecular basis of DLM phenotypes and to contribute knowledge about mechanisms of host-pathogen interaction.

New Insights on Rotenone Resistance of Complex I Induced by the m.11778G>A/MT-ND4 Mutation Associated with Leber's Hereditary Optic Neuropathy.

The finding that the most common mitochondrial DNA mutation m.11778G>A/MT-ND4 (p.R340H) associated with Leber's hereditary optic neuropathy (LHON) induces rotenone resistance has produced a long-standing debate, because it contrasts structural evidence showing that the ND4 subunit is far away from the quinone-reaction site in complex I, where rotenone acts. However, recent cryo-electron microscopy data revealed that rotenone also binds to the ND4 subunit. We investigated the possible structural modifications induced by the LHON mutation and found that its amino acid replacement would disrupt a possible hydrogen bond between native R340 and Q139 in ND4, thereby destabilizing rotenone binding. Our analysis thus explains rotenone resistance in LHON patients as a biochemical signature of its pathogenic effect on complex I.

Modelling eNvironment for Isoforms (MoNvIso): A general platform to predict structural determinants of protein isoforms in genetic diseases.

The seamless integration of human disease-related mutation data into protein structures is an essential component of any attempt to correctly assess the impact of the mutation. The key step preliminary to any structural modelling is the identification of the isoforms onto which mutations should be mapped due to there being several functionally different protein isoforms from the same gene. To handle large sets of data coming from omics techniques, this challenging task needs to be automatized. Here we present the MoNvIso (Modelling eNvironment for Isoforms) code, which identifies the most useful isoform for computational modelling, balancing the coverage of mutations of interest and the availability of templates to build a structural model of both the wild-type isoform and the related variants.

Structure, dynamics, and function of SrnR, a transcription factor for nickel-dependent gene expression.

Streptomyces griseus, a bacterium producing antibacterial drugs and featuring possible application in phytoremediation, expresses two metal-dependent superoxide dismutase (SOD) enzymes, containing either Fe(II) or Ni(II) in their active site. In particular, the alternative expression of the two proteins occurs in a metal-dependent mode, with the Fe(II)-enzyme gene (sodF) repressed at high intracellular Ni(II) concentrations by a two-component system (TCS). This complex involves two proteins, namely SgSrnR and SgSrnQ, which represent the transcriptional regulator and the Ni(II) sensor of the system, respectively. SgSrnR belongs to the ArsR/SmtB family of metal-dependent transcription factors; in the apo-form and in the absence of SgSrnQ, it can bind the DNA operator of sodF, upregulating gene transcription. According to a recently proposed hypothesis, Ni(II) binding to SgSrnQ would promote its interaction with SgSrnR, causing the release of the complex from DNA and the consequent downregulation of the sodF expression. SgSrnQ is predicted to be highly disordered, thus the understanding, at the molecular level, of how the SgSrnR/SgSrnQ TCS specifically responds to Ni(II) requires the knowledge of the structural, dynamic, and functional features of SgSrnR. These were investigated synergistically in this work using X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, atomistic molecular dynamics calculations, isothermal titration calorimetry, and in silico molecular docking. The results reveal that the homodimeric apo-SgSrnR binds to its operator in a two-step process that involves the more rigid globular portion of the protein and leaves its largely disordered regions available to possibly interact with the disordered SgSrnQ in a Ni-dependent process.

Nickel import and export in the human pathogen Helicobacter pylori, perspectives from molecular modelling.

The uptake of essential metal ions and the ability to extrude them when their excess causes toxicity are crucial processes for all living beings. Nickel is a virulence factor for several human pathogens and in particular for the human gastric pathogen Helicobacter pylori because of its crucial role in the catalytic activity of two Ni-dependent enzymes, urease and hydrogenase. H. pylori requires efficient uptake mechanisms to import Ni(II) because of its scarcity in the human body, but the molecular details of Ni(II) homeostasis are not fully known. Here we offer a structural framework for the machinery of Ni(II) import/export in H. pylori, obtained through comparative modelling and macromolecular docking. The model structures reported in this perspective are initial steps towards the understanding of these processes at the molecular level and in the direction to exploit them to eradicate infections caused by this family of pathogens. The differences between the structural models obtained by using both the recently released neural network-based approach implemented in AlphaFold2 and a more classical user-driven modelling procedure are also discussed.

Definition of the Binding Architecture to a Target Promoter of HP1043, the Essential Master Regulator of Helicobacter pylori.

HP1043 is an essential orphan response regulator of Helicobacter pylori orchestrating multiple crucial cellular processes. Classified as a member of the OmpR/PhoB family of two-component systems, HP1043 exhibits a highly degenerate receiver domain and evolved to function independently of phosphorylation. Here, we investigated the HP1043 binding mode to a target sequence in the hp1227 promoter (Php1227). Scanning mutagenesis of HP1043 DNA-binding domain and consensus sequence led to the identification of residues relevant for the interaction of the protein with a target DNA. These determinants were used as restraints to guide a data-driven protein-DNA docking. Results suggested that, differently from most other response regulators of the same family, HP1043 binds in a head-to-head conformation to the Php1227 target promoter. HP1043 interacts with DNA largely through charged residues and contacts with both major and minor grooves of the DNA are required for a stable binding. Computational alanine scanning on molecular dynamics trajectory was performed to corroborate our findings. Additionally, in vitro transcription assays confirmed that HP1043 positively stimulates the activity of RNA polymerase.

Probing the transport of Ni(II) ions through the internal tunnels of the Helicobacter pylori UreDFG multimeric protein complex.

The survival of several pathogenic bacteria, such as Helicobacter pylori (Hp), relies on the activity of the nickel-dependent enzyme urease. Nickel insertion into urease is mediated by a multimeric chaperone complex (HpUreDFG) that is responsible for the transport of Ni(II) from a conserved metal binding motif located in the UreG dimer (CPH motif) to the catalytic site of the enzyme. The X-ray structure of HpUreDFG revealed the presence of water-filled tunnels that were proposed as a route for Ni(II) translocation. Here, we probe the transport of Ni(II) through the internal tunnels of HpUreDFG, from the CPH motif to the external surface of the complex, using microsecond-long enhanced molecular dynamics simulations. The results suggest a "bucket-brigade" mechanism whereby Ni(II) can be transported through a series of stations found along these internal pathways.

A Blueprint for High Affinity SARS-CoV-2 Mpro Inhibitors from Activity-Based Compound Library Screening Guided by Analysis of Protein Dynamics.

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein's druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.